Can Quantum Lattice Fluctuations Destroy the Peierls Broken 

Symmetry Ground State? 

William Barford^ and Robert J. BursilP 
^Department of Physics and Astronomy, 
University of Sheffield, Sheffield, S3 7RH, United Kingdom 
"^School of Physics, University of New South Wales, 
Sydney, New South Wales 2052, Australia 

Abstract 

The study of bond alternation in one-dimensional electronic systems has had a long history. 
Theoretical work in the 1930s predicted the absence of bond alternation in the limit of infinitely 
long conjugated polymers; a result later contradicted by experimental investigations. When this 
issue was re-examined in the 1950s it was shown in the adiabatic limit that bond alternation occurs 
for any value of electron-phonon coupling. The question of whether this conclusion remains valid 
for quantized nuclear degrees of freedom was first addressed in the 1980s. Since then a series of 
numerical calculations on models with gapped, dispersionless phonons have suggested that bond 
alternation is destroyed by quantum fluctuations below a critical value of electron-phonon coupling. 
In this work we study a more realistic model with gapless, dispersive phonons. By solving this 
model with the DMRG method we show that bond alternation remains robust for any value of 
electron-phonon coupling. 

PACS numbers: 71.30.+h, 71.10.-w, 71.38.-k 
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The effect of quantizing the nuclear degrees of freedom on the bond alternation was 
studied by Fradkin and Hirsch^ on models of noninteracting electrons. They used both 
renormalization group arguments and Monte Carlo simulations in their investigations. In a 
model with gapless, dispersive (or Debye) phonons they showed that for any value of the 
electron-phonon interaction the metallic state is unstable with respect to a lower symmetry 
insulating phase that exhibits bond alternation. This result confirms the adiabatic predic- 
tion. The transition from the metallic state to the insulating state is called the Peierls 
transition, with the ground state named the Peierls state. They also predicted that for 
the same model with spinless fermions the Peierls state is destroyed by quantum lattice 
fluctuations below a critical value of the electron-phonon interaction. Since a model of spin- 
less fermions is equivalent to the XY quantum antiferromagnet on open chains, this result 
suggests that strong electronic interactions coupled to quantized phonons might destroy the 
Peierls state. 

More recent numerical calculations on various models have also indicated that the Peierls 
state may be destroyed by quantized phonons below a critical value of the electron-phonon 
interaction. Using a finite size scaling analysis of the spin gap Caron and Moukouri^i showed 
that at a critical value of the electron-phonon interaction there is a Kosterlitz-Thouless 
transition in the XY spin-Peierls model with gapped, non- dispersive (or Einstein) phonons. 
Similar conclusions were made for the Heisenberg spin-Peierls model, again with Einstein 
phonona^ii. The Holstein model with Einstein phonons for both spinlesa^ and spinfull^ 
fermions have also be shown to exhibit a Peierls phase transition at a non-zero value of the 
electron-phonon interaction. 

Of more direct relevance to the present work is the investigation by Sengupta et al.^ 
In this work the electronic degrees of freedom are modelled by a tight-binding model with 
onsite and nearest neighbour Coulomb interactions (namely, the extended Hubbard model). 
The nuclear degrees of freedom are described by Einstein bonds phonons. As in previous 
work with Einstein phonons, there is a critical value of the electron-phonon interaction below 
which lattice fluctuations destroy the Peierls state. 

In this article we focus our attention on the role of quantized gapless, dispersive (or 
Debye) phonons on the Peierls state. The electronic degrees of freedom are modelled by the 
extended Hubbard model, with the electronic and nuclear degrees of freedom being coupled 
via linear distortions of the bond lengths. We study this model as function of the electron- 
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phonon interaction, the Coulomb interaction and the phonon frequency. In all cases we 
consider a half filled band. In the limit of large Coulomb interactions this model at half- 
filling maps onto the quantum Heisenberg antiferromagnet, thus establishing a connection 
with our previous work®. Alternatively, in the limit of vanishing phonon frequency (namely, 
M oo) the model maps onto the classical, adiabatic model, which has been extensively 
studied (see, for example^ and references therein). 

The Hubbard-Peierls model with dispersive, gapless phonons is defined by, 

H = Hi + H2 (1) 

where 

Hi= - X] (i + Oi{ue - ue+i)) (cJ^Q+i^ + 
e=i,a 

+uf2(Nn-l) (iV,i-i) +i5^m.-l)(iVm-l)- (2) 



= + (3) 



and 

2M 2 

Ui and Pe are the conjugate displacement and momentum operators for the ith site, c^^ 
creates an electron with spin a on site i. t is the electronic hybridization integral, U and V 
are the onsite and nearest neighbour Coulomb interactions, respectively, a is the electron- 
phonon coupling constant, M is the nuclear mass, and K is the elastic spring constant. 

In the adiabatic limit (defined by M oo) and the noninteracting limit (defined by 
U = V = 0) the bond alternation, defined by 



N 

t 

(where a is the average bond length), satisfies 

(5o = 4 exp 



(4) 
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for A <^ 1. A is the usual definition of the electron-phonon coupling parameter, defined as 
A = 2a'^/7rKt. A non-zero value of the electron-phonon interaction. A, therefore implies a 
non-zero value for the bond alternation. The purpose of this work is to address the role of 
quantized phonons on the bond alternation, and in particular to investigate whether there 



is a non-zero value of A for the Peierls transition. We find that although quantized phonons 
reduce the amplitude of the broken symmetry order parameter, in contrast to the previous 
investigations using Einstein phonons, there is no evidence that the Peierls state is destroyed 
at a finite value of the electron-phonon interaction. 

I. METHODS 

The lattice degrees of freedom are quantized for dispersive phonons by introducing the 
phonon creation and annihilation operators, 



Eq. is solved by the density matrix renormalization group (DMRG) method^*i2iiiJ^ 
for chains of up to 136 sites. The DMRG method is an efficient truncation procedure 
for solving quantum lattice Hamiltonians, especially in one-dimension^^,. Solving electron- 
phonon models, however, poses special problems as the number of phonons is not conserved. 
Thus it is necessary to retain enough oscillator levels per site (Mp) to ensure convergence 
of the key quantities. The number of electron-phonon states per site (4 x Mp) is usually 
too many to augment with the system block, thus it is necessary to truncate this single site 
basis. There are therefore four convergence parameters to tune for this model: the number 
of oscillator levels per site, the number of optimized electron-phonon states per site, the 
number of states per system block, and the overall number of superblock states. 

We first establish convergence with respect to the number of optimized electron-phonon 
states per site. Table HI shows that both the ground state energy and the phonon order 
parameter (defined by Eq. Q) have converged to better than 0.2% for 14 states per site. 
Next, we consider convergence with respect to the number of system and superblock states. 
Table HTl shows excellent convergence when the number of system block states is 250 and the 
eigenvalue product cutoff e = 10~^^. Finally, we consider convergence with respect to the 
number of bare oscillator levels per site. Table HTTl and Table HVl shows excellent convergence 
at 6 levels per site for both the ground state energy and the phonon order parameter. 




(6) 



and 




(7) 



respectively, where uj = ^J^KJM^^^. 
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TABLE I: The ground state energy, E, (in units of t) and the phonon order parameter, Nm, 
(defined by Eq. Q) of the Hubbard-Peierls model as a function of the number of electron-phonon 
states per site, M, for a 40-site chain with 6 oscillator levels per site, to = t, U = 2.5t, V = U /4, 
and g = 0.1. 

M E Nm 
6 -58.716 0.663 
10 -59.969 0.629 
14 -60.220 0.628 
18 -60.275 0.629 



TABLE II: The ground state energy, E, (in units of t) and the phonon order parameter, Nm, of 
the Hubbard-Peierls model as a function of the density matrix eigenvalue product cutoff, e, the 
number of system block states, Af, and the superblock Hilbert space size, SBHSS, for a 40-site 
chain with 6 oscillator levels per site, uj = t, U = 2.5t, V = J7/4, and g = 0.1. 
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TABLE III: The ground state energy (in units of t) of the Hubbard-Peierls model as a function of 
the number of sites, A^. to = t, U = 2.5t, V = U /4, and g = 0.1. 
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TABLE IV: The phonon order parameter, Nm, of the Hubbard-Peierls model as a function of the 
number of sites, N. uj = t, U = 2.5t, V = U /A, and g = 0.1. 

Number of oscillator levels per site 

N I 2 3 4 5 6 7 

16 0.321 0.367 0.367 0.384 0.385 0.385 

24 0.398 0.455 0.455 0.477 0.479 0.479 

40 0.518 0.593 0.595 0.625 0.628 0.629 



Typically, we use 6 oscillator levels per site, with 14 optimized single site states, 250 
block states, and 500, 000 superblock states. One finite lattice sweep is performed at the 
target chain size. In addition, in situ optimization of the single site basis is performed for 
each site at the target chain^- . In all cases we investigate linear chains with open boundary 
conditions. We maintain constant chains lengths by fixing the position of the end sites. 

The Hubbard-Peierls model in the classical limit (defined by Eq. is solved using the 
Hellman-Feynamn theorem which implies that, 

- ut) = (^{cla^e+ia + 4+1^^^) - , (8) 

where (■ ■ ■ ) means the expectation value in the ground state. The ground state bond 
alternation is found by iterating Eq. (jH)) until there is convergence^»i^. 



II. RESULTS AND DISCUSSION 



We first investigate the Hubbard-Peierls model for arbitrary electron-phonon interactions 
in the noninteracting limit {U = V = 0). Following Fradkin and Hirsch^, we use the 
dimensionless electron-phonon interaction, g = a/ / Kt = a/Avt/^. 

The Peierls broken-symmetry ground state is characterized by a non-zero value of the 
staggered phonon order parameter, defined by 

where Bi = {b\ + be)/2 is the dimensionless displacement of the ith site and is the number 
of sites. Notice from Eq. (jH)) and Eq. ((Tj) that the bond alternation, 60, is proportional to 
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FIG. 1: The phonon order parameter versus the electron-phonon interaction for the noninteracting 
Hubbard-Peierls model, uj = t. The inset shows a Hnear extrapolation of Nm versus g to the 
origin. (The crosses indicate the evaluated points.) 

^/uJm{N). To determine whether this order parameter vanishes as a function of the electron- 
phonon interaction in the asymptotic limit it is necessary to perform a finite size scaling 
analysis. Fradkin and Hirscb^ suggested that the scaling relation 

m{N) = -^F{Nm{oo)) (10) 

should apply. Thus, curves of Nm{N) versus g are expected to coincide at the critical value 
of g at which m(oo) vanishes. 

Fig. (0) shows Nm versus g for values of uj = t and U = V = 0. These results were ob- 
tained using the density matrix renormalization group method, as described in the Methods 
section. It is clear that for small values of g m is linearly proportional to g for all chain 
lengths. Furthermore, the curves extrapolate linearly to the origin indicating a zero value for 
the critical g. This confirms Fradkin and Hirsch's original prediction for the noninteracting 
limit of this model-, and accords with the adiabatic prediction. Notice that since the distor- 
tion of the ith bond from its uniform value, (m£+i — M£+i), is proportional to ^JuJg{BiJ^l~B(), 
the bond alternation is oc or oc A for g <^1. 

We next consider the role of Coulomb interactions on the value of the bond alternation. It 
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FIG. 2: The phonon order parameter versus the electron-phonon interaction for the interacting 
Hubbard- Peierls model, uj = t, U = 2.5t, V = U /A. The inset shows a hnear extrapolation of Nm 
versus g to the origin. (The crosses indicate the evaluated points.) 

is known that in the adiabatic limit the bond alternation initially increases as the Coulomb 
repulsion increases. This is because electronic interactions suppress the quantum fluctuations 
between the two degenerate bond alternation phases^J^. The bond alternation is maximized 
when the electronic kinetic energy roughly equals its potential energy, namely when U ~ 4t. 
As the Coulomb interaction increases further the bond alternation decreasesS*^, eventually 
scaling as t/{U — V). 

Fig. shows Nm versus g for values of = t, [/ = 2.5t, and V = U/4. These Coulomb 
parameters were chosen to coincide with those used by Sengupta et al^ , who studied a similar 
model to Eq. ((T)), but with Einstein bond phonons. As already stated, they reported that 
quantized phonons destroy the Peierls ground state below a critical value of the electron- 
phonon interaction. In contrast to that work^, however, we see that for our model with 
Debye phonons the predictions are qualitatively similar to the noninteracting limit, namely 
Nm decreases as a linear function of g a.s g —>■ 0. We note, however, that the phonon order 
parameter for a given value of g is larger for this value of U than for the noninteracting 
limit. 
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FIG. 3: The phonon order parameter for 72-site chains versus the Coulomb interaction U for the 
Hubbard- Peierls model, to = t, V = U/A, and g = 0.01. 

The phonon order parameter, m, as a function of the Coulomb interaction U is illustrated 
in Fig. (j^I) for 72-site chains, g = 0.01, and thus the results lie in the linear regime where m oc 
g. As for the adiabatic limit, m initially increases with U before asymptotically decreasing 
as U ^ oo. We might therefore expect that in the strong coupling limit {U oo) quantum 
fluctuations have a better chance to destroy the Peierls state below a critical electron- 
phonon interaction. That this is not the case was demonstrated by the present authors in 
an investigation of the Heisenberg spin-Peierls model^. This model is the Hubbard-Peierls 
model in the strong coupling limit at half-filling. The Heisenberg spin-Peierls model is 
parameterized by J = 4t^/([/ — V),uj and g. Re£^ showed that as for the Hubbard-Peierls 
model, Nm also decreases as a linear function of g a.s g ^ 0. We have therefore shown that 
for all electronic interactions the Peierls state is robust to quantized lattice fluctuations. 

We now turn to discuss the phonon order parameter as a function of the phonon frequency, 
uj. As uj (and more oscillator levels are occupied) the system approaches the classical 
limit. To make a direct comparison to the classical limit we also solve Eq. in this limit. 
The classical Hamiltonian is defined by. 



(where Hi is defined by Eq. and the bond variables are now classical variables. The 




(11) 
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FIG. 4: The scaled bond alternation of the Hubbard-Peierls model for a 50-site chain. The adiabatic 
limit (using Eq. ((TT|)) (diamonds); the quantum case (using Eq. w = O.lt (circles), and oj = t 
(triangles). U = 2.5t and V = U/A. 

role of the last term is to ensure constant chain lengths, where F is self-consistently chosen 
so that ^^(-Uf — M^+i) = 0. The ground state bond alternation of the classical model is found 
iteratively using the Hellmann-Feynman procedure, as described in the Methods section. 

Fig. (jU shows 6q (scaled by g) of a 50-site chain versus g for the two quantum cases, 
u = O.lt and u = t, as well as for the adiabatic limit. Evidently S^/g tends linearly — * as 
g ^ 0. Thus, 6o Qc g^ Qc \. Moreover, as predicted, the phonon order parameter decreases as 
a function of cj. As g increases the deviation between the quantum and classical predictions 
increases. 

Finally, we discuss the nature of the Peierls transition at A = 0. In the noninteracting 
limit the classical transition is given by Eq. (jSJ. No analytical expression exists in the 
interacting limit. However, Fig. Q shows that the transition (at g = 0) is qualitatively 
similar in the classical and quantum cases, as 6o is proportional to g^ = A. Moreover, 
the behaviour of the bond alternation as a function of chain length is the same for small 
electron-phonon interactions in the classical and quantum cases. Fig. (0) shows the scaled 
bond alternation versus inverse chain length in the adiabatic and quantum cases. The close 
agreement of these predictions suggests that the transition at (7 = is the same in the 
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FIG. 5: The scaled bond alternation versus inverse chain length in the adiabatic limit. Crosses are 
the quantum result with uj = t. g = 0.1, U = 2.5t, and V = U/4:. 

quantum and classical cases in the limit of long chains. 

III. CONCLUSIONS 

In conclusion, we have used the DMRG method to investigate the role of gapless, disper- 
sive quantized phonons on the Peierls transition in the Hubbard-Peierls model. We showed 
that the staggered phonon order scales as g (and the bond alternation scales as g"^, and thus 
oc A) as (7 — > for all values of the Coulomb interaction and phonon frequency considered 
here. In the strong coupling limit the Hubbard-Peierls model maps onto the Heisenberg spin- 
Peierls model, which also exhibits a Peierls transition at g = 0^. Thus, we conclude that 
the Peierls transition occurs at g = for all values of the Coulomb interaction. Moreover, 
we showed that the quantum predictions for the bond order follow the classical prediction 
as a function of for small g. We therefore conclude that the zero g phase transition is 
of the mean-field type. 
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